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An f(R) gravity model is proposed to realize a late time accelerated expansion of our Universe. To test the 
viability of an f(R) gravity model through cosmic observations, the background evolution and the Einstein- 
Boltzmann equation should be solved for studying the effects on the cosmic microwave background power 
spectrum and on the matter power spectrum. In the market, we already have the modified versions of CAMB 
code, for instance EFTCAMB and MGCAMB. However, in these publicly available Einstein-Boltzmann codes, 
a specific background cosmology, for example the ACDM or wCDM, is assumed. This assumption would be 
non-proper for a specific f(R) model where the background evolution may be different from a ACDM cosmol¬ 
ogy. Therefore the main task for this paper is to present a code to calculate the anisotropies in the microwave 
background for any f(R) gravity model based on CAMB code, i.e. FRCAMB, where the background and per¬ 
turbation evolutions are included consistently. As results, one can treat FRCAMB as a blackbox to output the 
CMB power spectrum and matter power spectrum, once an f{R) function, its first two derivative with respect to 
R, i.e. f R = df/dR, f rr = cPf/dR 2 and the reasonable values of the model parameters are inputted properly. As 
by-products, one can also output the effective equation of state of f(R ) model, the evolution of the dimensionless 
energy densities and other interesting cosmological quantities. 


I. Introduction 

The late time accelerated expansion of our Universe de¬ 
mands a modification of general relativity (GR) at large scale 
or an addition of an extra exotic energy component, see the 
monograph [1] and references therein. Usually one believes 
that a modified gravity (MG) model is degenerate to a dark 
energy model at the background level. That is to say for a 
modified gravity model there exists a dark energy model hav¬ 
ing an effective equation of state (EoS) which produces the 
same expansion history as that of the modified gravity model. 
Thus for the same background evolution one expects to distin¬ 
guish modified gravity from dark energy models through the 
dynamical evolutions of perturbations. In the literature, many 
efforts have been made aiming to detect the possible deviation 
to GR through a parameterized Poisson equation and the slip 
of the Newtonian potentials while fixing the background evo¬ 
lution to a ACDM cosmology. The parameterized formalism 
with a fixing ACDM background is proper, when one only 
concerns a general MG model which probably has no explicit 
Lagrangian in general. Thus the fact would be embarrassed 
even if a significant deviation from GR is confirmed, because 
the explicit MG theory is still unknown in this parameterized 
formalism. Of course before fixing a proper MG, a significant 
deviation to GR should be confirmed by cosmic observations. 
This is the main reason why this parameterized formalism is 
still plausible now. However unfortunately currently available 
cosmic observations have not found any significant deviation 
to GR [2], 

In this paper, instead of considering the parameterized for¬ 
malism, we will mainly focus on a family of modified gravity, 
i.e. f(R) gravity model, not only including the background 
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evolution but also including the perturbation evolution. When 
an f(R) gravity model is considered, it becomes non-proper to 
keep a ACDM background evolution. It is mainly because that 
a ACDM background evolution specifies the form of f(R) [3]. 
As well known, a reasonable Universe should mainly experi¬ 
ence three important epochs consequently: radiation, matter 
and (effective) dark energy dominated stages. In terms of the 
effective EoS 


our Universe should vary from w e ff = 1/3 (radiation domi¬ 
nated epoch), w e ff = 0 (matter dominated epoch) to w e ff < 
-1/3 (dark energy dominated epoch) at different epochs. 
Therefore the background evolution provides a preliminary 
and valuable test for the viability of an f(R) gravity model. 
The geometrical measurements from the luminosity of type 
la supernovae (SNe) as standard candle, the angular diame¬ 
ter distance of baryon acoustic oscillation (BAO) as standard 
ruler and the positions of the peaks of the cosmic microwave 
background (CMB) radiation power spectrum can constrain 
an f(R) gravity model at the background level extensively. 

In addition to the geometrical measurements, the dynamical 
measurements related to the evolution of the perturbation at 
the linear and non-linear scales break the degeneracy of model 
parameters and provide even stronger constraint to an f{R) 
gravity model. To study the CMB power spectrum, one should 
study the perturbation evolution for an f(R) gravity model, 
i.e. solve the full Einstein-Boltzmann equation for photons. 
In the market, we already have the Einstein-Boltzmann equa¬ 
tion solvers for the modified gravity models (including f(R) 
gravity model as a special case), for instances MGCAMB 
[4], EFTCAMB [5] and FRCAMB [3] which are modified 
version of CAMB package [6]. However for the MGCAMB 
code, a ACDM or wCDM background is assumed when the 
CMB power spectrum is shown in an f(R) gravity model. Ob¬ 
viously it is non-proper once the background evolution for an 
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f(R) gravity model deviates from that of ACDM or wCDM. 
The situation is little different for the EFTCAMB code, where 
a background evolution is fixed and an f(R ) model is recon¬ 
structed from this fixed expansion history. Although the back¬ 
ground evolution is consistent to the perturbation evolution, 
the freedom of choosing different f(R) gravity model is lost. 
Our FRCAMB code is designed for any f(R) gravity model at 
the background and perturbation evolutions based on the pub¬ 
licly available FRCAMB code [3], once a form of f(R) and 
its first two derivatives with respect to R are inputted. 

This paper is structured as follows. In Section II, the back¬ 
ground evolution equations for an f{R) gravity model are pre¬ 
sented. The perturbation evolution is given in Section III. The 
CMB power spectrum and matter power spectrum are given 
in Section IV. Section VI is the conclusion. 


and pi, i = m, r is the energy density of the matter (cold dark 
matter+baryon) and radiation. Replacing the Ricci scalar R 
by R - 12 H 2 + 6 HH', the above equation (5) can be recast to 

H 2 - f R (HH' + H 2 ) + f - + H 2 f m R' = mV 3 + ^a~ 4 ), (6) 
6 Li„, 

where 

m 2 = “ (8315 Mpc)- 2 , (7) 

and Q, = K 1 p l o/?>II 2 , i = m, r is the dimensionless energy den¬ 
sity of the matter (cold dark matter+baryon) and radiation. 
The subscript 0 denotes the corresponding value at present 
a = 1. In similar to [11], defining the dimensionless variables 


II. Background Evolution for an f(R) Gravity Model 

The Einstein-Hilbert action for an f(R) gravity model reads 
as 




H 2 

m 2 

R 


n r 


—4 


y R = — - 3a 


( 8 ) 

(9) 


5 f d4x ^[R + /(*)] + f d 4 x y/—g(-C m + £ r ), 

( 2 ) 

where and -£, are the Lagrangian of matter and radiation 
respectively, which will not include the mysterious dark en¬ 
ergy as the late time accelerated expansion of our Universe 
can be realized by the proposed f(R ) gravity. For recent re¬ 
views for modified gravity theory, see Refs. [7-10]. Doing 
variation with respect to the metric g tIV for the Einstein-Hilbert 
action, one obtains a generalized Einstein equation which re¬ 
lates the geometry of space-time to the distribution of energy- 
momentum 

+ f R R MV - iy-f - □/*) gtlv - V„V v /« = 8 t tGT^, (3) 


and using Eq. (6), one obtains the differential equations for 

{yH,y R } 


Vh = if* - 4 fH, 


( 10 ) 


y' R = 9a 3 - 


m 2 f RR (vh + a 3 + 
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r if* 1 -3 

}’H - f R \ — - y’H - -a 


a _ 4 

—a 

1 in. 


6m 2 


.( 11 ) 


The solutions of this equations describe the background evo¬ 
lution for an f(R ) gravity model, once the initial conditions at 
present ao = 1 are given 


where 


/* = 


df(R) 

dR 


(4) 


plays as an effective scalar filed. It is obvious that the gen¬ 
eral relativity is recovered when f(R) - 0. The Friedmann 
equation for an f(R) gravity reads as 


H 2 + f R (h 2 - + f - + H 2 f RR R' = j {p m + Pr ), (5) 

where k 2 = d>nG, f RR = d 2 f{R) / dR 2 , the prime ' denotes the 
derivative with respect to the nature logarithm of the scale fac¬ 
tor a, i.e. In a, H — a/a is the expansion rate of our Universe; 


y^ a = 1 m 2 1 q ’ 

III 

H l 

f*Ui = 6(1 - go)—I - 3, 


( 12 ) 

(13) 


where we have used the relation Rq = 6// 3 ( 1 - qo), here qo is 
the value of the deceleration parameter q = -aa/cr at present. 
<7o is a free model parameter which is related to f R o when an 
f(R) gravity model is known. Therefore as a comparison to 
the literature, f R o is a derived model parameter. 

The effective dark energy density and pressure for an f{R) 
gravity are given by [1] 
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k 2 Pde = 2 VrR - /) - 3H 2 f RR R' + 3 // 2 (/*> - /r) 


" 2 ^ (/r4 - 4) - 3^/rr*' + 3^ (/ro - /r) 


K PDE 


= f RR H 2 + /**/?'//#' + 2H 1 f R - I (*/* - /) - (2 HH' + 3H 2 ) (/*> - /r) 


-r 

2 VR 

iii- 


1 iH 1 


2 \ m 2 


. //- 


- — +2 — 


1 




/ 


/RR*' - - /r—-T I — 


H 1 


--A + 3^ 

m z / ra z 


(/ro ~ /r) f ■ 


(14) 


(15) 
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The EoS of the effective dark energy wqe = Pde/Pde can be 
obtained easily 


Woe - -1+ 


2 {/rH 2 + f R H'H) - 2 H 2 f' R - 4 HH' (f m - f R ) 


= -1 


/r* - f 

J l\ ,11 - 111 - 


6 Kf RR R' + eK (jRo-f R ) 


The effective EoS of our Universe becomes 


2 h ' i(gy 

Weff 3 H 3 (fR) ’ 


(17) 


which enables us to test the viability of the f(R) model at the 
background level. The main reason is that w e ff should vary 
from w e ff = 1/3 (radiation dominated epoch), w e ff = 0 (mat¬ 
ter dominated epoch) to w e ff < -1/3 (dark energy dominated 
epoch) at different epochs. With the help of the following ex¬ 
pressions, 




f RR R' = m 2 f RR (y R - 9a 3 ), 
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derivatives of Eq. (22) with respect to R easily 
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(23) 


(24) 


Only these three functions f(R), f R and f RR are needed in our 
FRCAMB package. In order to make the expansion history 
close to that of ACDM [11], it was shown that the parameters 
(ci, C 2 ) were related to £ 2 „,o, O,o 


Cl (1 - 1 r{) O m o) 

— = o- 


c 2 


£1 


(25) 


m() 


H 1 


1 


Q, 


— I = 3 ^ - 4 y« - 3a - 4 —a 


finally the Eqs. (14), (15), (16) and (17) can be rewritten as 
functions of the dimensionless variables { yH,yR )■ Until now, 
we have obtained the background evolution 


H 2 = 


3(1+ /ro) 


(Pm + Pr + PDE) ■ 


Therefore for the HS model, we have three free model param¬ 
eters {ci, qo, n). In this work, we consider the n = 1 case only. 
It is easy to extend to the n + 1 cases. We show the relative 
difference Uf flR) (z) - H GR (z))/H GR (z ) x 100% for the HS and 
ACDM model in Figure 1 by adopting the same cosmolog¬ 
ical parameters obtained by Planck [12] but with varying ci 
or effectively varying log(|/ M ) - 1|). The results show that the 
background evolution of HS model for n - 1 is very close to 
+ 4.2lL a ( fjjj /hat of the ACDM. The relative deviation to the ACDM back- 
£/» ground is less than 0.07%. This also confirms the viability of 
( 19 pur code for the background evolution. The curves for dif¬ 
ferent values of c\ imply that the background evolution is not 
( 2 ()fcnsitive to the values of c\ for the HS model in n — 1 case. 
We also show the evolution of the dimensionless density pa¬ 
rameter Clx=m,r,de(z) and the effective EoS with respect to the 
redshift z in Figure 2, where one can see a series of transitions 
of our Universe from the early radiation dominated epoch, the 
middle dark matter dominated epoch to the late effective dark 
energy dominated epoch. It says that the HS model pass the 
background evolution test. 


( 21 ) 


for an f(R) gravity and the effective energy density, pressure 
and EoS. 

Taking the Hu-Sawicki (HS) model 


m = m 2 


R 


c 1 


( 5 )" 




( 22 ) 


Taking the trace of Eq. (3), one has 

3□/*-/? + fuR - 2 / = -K 2 (p m + Pr). 


(26) 


This equation can be recast as a evolution equation of a scalar 
field f R 


□/r 


Weff 
df R ’ 


as a working example in this paper, one obtains its first two 


(27) 
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FIG. 1. The relative difference, (Hf^(z) - H GR {z))/H aR (z) x 100%, 
of the Hubble parameters for the HS (n = \.qa = -0.65) and ACDM 
model. 


Usually, instead of using the Compton wavelength Af R = , 

one defines the dimensionless Compton wavelength 


JRrR' H _ fRRR' „,2 
T7TrW~ 1+fRUHiV- 

2 V m 2 ) 


Taking the HS model (n = 1) as a working example, we 
show the evolution of /«, m 2 R and B with respect to the red- 
shift z in Figure 3, where m 1 ^ begins to increase quickly at a 
transition point dependent on the values of model parameter. 
It implies that the f(R ) gravity approaches to GR when the 
effective scalar field /r becomes massive. Therefore to make 
the code efficient and stable, we put a cutoff to the mass of 
m 2 ^ ~ 5 x 10 2 Mpc 2 at the scale factor a - a, which value 
is determined by the model parameters. For the small val¬ 
ues of a < a t , we replace the evolution and perturbation with 
the standard GR model with the same cosmological model pa¬ 
rameters. Actually in the early times a < a,, the contribution 
of the effective dark energy component can be neglected as 
shown in Figure 2. 



FIG. 2. The evolution of the dimensionless density parameters of 
flx=m,r,de(z) an d the effective EoS with respect to the redshift z for the 
HS (n = 1. q 0 = -0.65 and log(|/«„ - 1|) = -5.75). 


with the effective potential 

dVgff 1 r 2 t 

= -\R-f K R + 2/ - K 2 (p m +p r )\. (28) 

which has an extremum at 



FIG. 3. The evolution of f R (blue curves), m 2 R (red curves) and B 
(green curves) with respect to the redshift z for the Hu-Sawicki model 
n = 1 with different values of ci or the effective log(|/Ro - 1|) = 
-5.75, -4.75, -3.75, -2.75. m 2 R begins to increase quickly at a 
transition point a, dependent on the values of ci or the effective 
logO/k) - 1|) = -5.75,-4.75,-3.75,-2.75. 


R-f R R + 2f = K 2 (p m +p r ). (29) 


The effective mass of this scalar field at the extremum is given 
by 



d2 Veff 

df 2 R 


1 

3 


1 +/* 

/«« 



(30) 


III. Perturbation Equations for an f(R) Gravity Model 


In this section, we will consider the scalar and tensor pertur¬ 
bation equations for an f(R) gravity which has already been 
well studied in Ref. [13, 14]. The line element with the scalar 
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and tensor perturbation can be written as 


and 


ds 2 = a 2 [-( 1 + 2 x VY (s) )dT Z + IBYfdTdx' 

+ (1 + 20 Y is) )j ij dx i dx j + SY < i j > dx‘dx j 
+ 2h T Y ( P dx' dx j ], (32) 

where y, ; - the three-dimensional spatial metric in the spherical 
coordinate is written as 


lYiji 


1 -Kr 2 

0 


0 0 

r 0 

0 r 1 sin 2 0, 


(33) 


and T w , Y ( f\ Y ( f and yP are the scalar and tensor harmonic 

J ’J 'J 

functions defined by 


(A + k 2 )Y (s) = 0, 


yC 5 ) —_y(s) 


Y K 


is) _ 


k ’ 

/ 0 s ) 


=_y (5) + _ Y ..y(*) 

“a 2 10 3 r,J ’ 


(34) 


(A + k 2 )YP = 0. 


In the synchronous gauge, by setting O' = 0 and II - 0 and 


k 2 

6F" + 2<H.5F' + a 2 (-r + ml )6F, 
a - JR 

= ^(Sp-36p)- ] -F’h’. (43) 

In this section, we have used the notation F = 1 + f R and 
the superscript' = d/dr. 'll = a! fa is the conformal Hubble 
parameter. 

In the CAMB package, the curvature perturbations are 
characterized by Z and <x 



cr = ka. 


where 


rf T = \<L<r-Z). ( 44 ) 

With the above variables, the perturbed modified Einstein 
equations recast into 


& 

q T = -(® + -), 

o 

(35) 

II 

On 

e 

(36) 


where qj refers to the conformal 3-space curvature perturba¬ 
tion 


6R m = 6 6K = -4(k~ - 3K)jj T , 


(37) 


the perturbed modified Einstein equations in the synchronous 
can be written as 


-ica 2 6p = -{-FFi + -F')h', - -'HSF’ - -6Fk 2 
2 2 4 L 2 2 


+ Fq T (k 2 - 3 K) + -FC5F , 

K 2 a 2 5p = F[-^Hh' L + | k 2 T) T - \h'[ - 2 ti t K] 


(38) 


+ 6F[F( 2 


a 


a 3 
5F" - 6F'Ft, 


k 2 + 2K] - \f'K l 


, VU O. F ' 2 2 P n SF 
a = -2Ma + m -ya-Ka - —, 


k 2 -3K , 1 , , 1 , 1 , F/j /f 

- Fq' = —K 2 a 2 q + -MF - -IYH5F + —. 

k T 2 2 2 2k 


(39) 

(40) 

(41) 


where 


_ (h L + 677r)' 

“ “ 2k 2 ’ 

q - (p + p)v, 


( 42 ) 


K 2 3 

(FV. + -F')kZ = — a 2 bp + Fk 2 q T /3 2 ~ -V.5F' 

- -6Fk 2 + -Ft’dF, (45) 

2 2 

y F(y6 2 o- - Z) = y crq + ^k6F ' 

- -k'HSF, (46) 

2 


cr'+2'Hcr-\ cr-kqj-K 2 a 2 - -k—, (47) 

F ' FkF 

z+( ^ + ^ )Z={ ^ + i + ^z 

2 F y 2 k F 

tea 2 3 5F" 

__(, p + 3 ^)--_ (48) 


where 


02 = 



(49) 


is the curvature factor. The propagation of the perturbed field 
<JF is given by 


k 2 

6F" + 2<H6F' + a 2 (— + ml )5F 
a 1 JR 

2 2 

= —y-(<5p - 36p) - kF'Z- (50) 

The source term of the CMB temperature anisotropy is 
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given by [15, 16] 


Sr(T,k) 

= e- £ (a" + rf T ) 

V, r 

+ g(A-po + 2a' H—4-—— t- —— 

k 12^57 4k 2 yffa 

+ g'(a 4—- + , . 

k 2k 2 yljh 4 k 2 


r 


r , , i g"c 
-) ■ 


_ E cr” kcr kZ, 
= e e (- + - -—) 


(51) 


cr' v r r" 

4- g(A T o + 2 — 4- — + ——— 4- ———) 

* k k nyjfc 4 k 2 ^_ 

, , ( <r Vb C ) | 1 g”( 

8K k k 2k 2 y/Jh 4 k 2 y/fc’ 

where g = —se~ B — an e (XTe~ e is the visibility function and e 
is the optical depth, f is given by 


3 9 

f = ( 4/2 + 2 ^ 2 ), 


(52) 


where I 2 , £2 indicate the quadrupole of the photon intensity 
and the E-like polarization respectively [16]. 

The propagation of gravitational waves for an f(R ) gravity 
is given by 


hj + 2'Tf 11 4- - 


1 dlnF\ 


2 din a ) 


h'+cUk 1 + 2K)h T = 


2>7iGcrW 


(53) 


where c 2 f is the square of the speed of gravitational waves. 

For the perturbation equations, one can clearly see that the 
only f(R ) gravity model dependent term is m 2 ^ in the per¬ 
turbed field equation (50) for 5F. For a concrete f(R) form, 
m can be calculated easily from its first two derivatives of 
f(R ) with respect to R. 


IV. CMB Power Spectrum and Matter Power Spectrum 

In this section, we will show the CMB power spectrum 
and the linear matter power spectrum for an f{R) gravity 
model. As mentioned in the previous section, the only in¬ 
puts for an f(R) gravity model are the three functions /(£), 
fn = df(R)/dR and /«« = d 2 f(R)/dR 1 . By setting the proper 
values of the model parameter for an f(R ) gravity model, 
our code FRCAMB will calculate the background evolution, 
solve the Einstein-Boltzmann equation and output the CMB 
power spectrum, the matter power spectrum, the effective EoS 
w e ff(a) and the dimensionless energy density ilx(a), and al¬ 
most all the quantities you are interested in. 

Taking the HS model (n — 1) as a working example, we 
show the CMB XX = TT,EE,TE,BB power spectrum in 
Figure 4, the effects on the CMB TT power spectrum and the 
linear matter power spectrum for different values of the model 
parameter ci or the effective logd /A) - 1|) in Figure 5 and Fig¬ 
ure 6 , where <70 = -0.65 is fixed for illustration. 

As shown in Figure 5, the CMB TT power spectrum is 
much sensitive to the values of the model parameter c 1 or its 


equivalent log(|/«o - 11 ) at low ( multipole where it is dom¬ 
inated by the late integrated Sachs-Wolfe (ISW) effect. This 
is mainly due to the fact that at the early epoch of our Uni¬ 
verse the mass of the effective scalar field /* becomes massive 
and there is no significant deviation to ACDM model, because 
ACDM cosmology is switched on when m 2 R ~ 5 x 10 2 Mpc ~ 2 
is arrived as shown in Figure 3. The same effects can be seen 
in Figure 6 , the linear matter power spectrum is sensitive to 
the values of c 1 or its equivalent logd/so — 1 |) after the matter- 
radiation equality epoch. 



FIG. 4. The CMB TT, EE, TE and BB power spectrum for HS model 
with n = 1, logd/ffo - 1|) = -5.75. 



FIG. 5. The ratio of the CMB TT power spectrum for different values 
of logd/flol) to that of the ACDM cosmology with the same values of 
the cosmological parameters. 
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FIG. 6. The linear matter power spectrum at the redshift z = 0 for 
different values of logO^o — 1|) and GR with the same values of the 
cosmological parameters. 


Once we have the linear matter power spectrum, the non¬ 
linear matter power spectrum can be calculated through a rou¬ 
tine like HALOFIT [17]. For this kind of halo-fit formula 
for an f(R) gravity model, an /V-body simulation is strongly 
demanded. We already have MGHalofit for calculating the 
nonlinear matter power spectrum [18], but it is only for the HS 
model (n = 1) and works in the range |/«)| e [ 10 6 , IGF 4 ] and 
z < 1 ■ A general halo-fit formula for any f{R) gravity is still 
unavailable. In our previous study on a specific family of f(R) 
gravity model, it is found out that the redshift space distor¬ 
tion (RSD) /cr§ can provide a tight constraint to the values of 
\/ro~ 1| ~ 1(T 6 [19]. For this small values of \fm— 1| ~ ICC 6 , 
the nonlinear matter power spectrum can almost mimic that 
of the ACDM model. With the very small values of /«o, it is 
difficult to detect a model not only because of the accuracy 
of the fitting formula but also the complicated astrophysical 
systematics on such scales [18]. 


V. Constraint to HS model from Cosmic Observations 

In this section, we show the constraint results to HS model 
for the n = 1 case from the geometric and dynamic measure¬ 
ments. For the geometrical one, we will use the supernova la 
data from SDSS-II/SNLS3 joint light-curve analysis [20], the 
baryon acoustic oscillation Dy(0. 106) = 456 + 27 [Mpc] from 
6 dF Galaxy Redshift Survey [21]; D v (035)/r s = 8.88 ± 0.17 
from SDSS DR7 data [22]; D v (0.51)/r s = 13.62 ± 0.22 
from BOSS DR9 data [23], the present Hubble parameter 
H 0 = 73.8 ± 2.4 [km s-'MpW 1 ] from HST [24], and the full 
information of CMB recently released by Planck2013 (which 
include the high-1 TT likelihood (CAMSpec) up to a maxi¬ 
mum multipole number of l max = 2500 from l - 50, the 
low-1 TT likelihood (lowl) up to / = 49) [25] with the addi¬ 
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1 ^ 7^9+0-037 
iJ,/ -0.03 7 
c ti +0.07 
' -0.16 
_1 qq+0.28 
a,: ^-0.63 

68.81 

0.7049 

0.2951 

o~ 8 

£re 


0.8168 

9.24 

Age/Gyr 

log(| fa - ID 
log(Bo) 


13.754 

-5.86 

-1.78 


TABLE I. The mean and best fit values with lcr errors for the inter¬ 
ested and derived cosmological parameters, where the Planck 2013, 
WMAP9, BAO, SN, HST and RSD data sets were used. 

tion of the low-1 TE, EE, BB likelihood up to / = 32 from 
WMAP9. For the dynamical one, we use the redshift space 
distortion (RSD) data. For using the growth rate, we calculate 
the fcr$(z) — dir?,/d In a at different redshifts in theory. For 
details, please see Ref. [19]. 

We perform a global fitting on the Computing Cluster for 
Cosmos by using the publicly available package CosmoMC 
[26] in the following model parameter space 

P = [£Vt 2 , Q c li 2 , 1000 mc, P n s , ln(10 10 A s ), q 0 , log(ci)}, 

(54) 

their priors are shown in the second column of Table I. The 
running was stopped when the Gelman & Rubin R- 1 param¬ 
eter /? — 1 — 0.02 was arrived; that guarantees the accurate 
confidence limits. The obtained results are shown in Table I 
and Figure 7. 


VI. Conclusion 

In this paper, we present an Einstein-Boltzmann equation 
solver, named FRCAMB, for calculating the anisotropies in 
the microwave background in any f(R) gravity model based 
on a modified version of CAMB. In this code, instead of as¬ 
suming a ACDM or uCDM background cosmology as done 
in the EFTCAMB code and the MGCAMB code, we solve 
the background evolution numerically for any f(R) gravity 
model once the f(R) function, its first two derivative with re¬ 
spect to R , i.e. /« s df/dR, f 'rr = d 2 f/dR 2 and the rea¬ 
sonable values of the f(R) model parameters are inputted. 
The outputs of this code include the CMB power spectrum, 
the matter power spectrum, the evolution of the total effec¬ 
tive EoS of our Universe and almost everything interesting. 
By global fitting to the model parameter space through the 
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Ml /so- 1 11 !) M B o) 


FIG. 7. The ID marginalized distribution and 2D contours for in¬ 
terested model parameters with 68% C.L., 95% C.L. by using the 
Planck 2013, WMAP9, BAO, BAO, JLA, HST and RSD data sets. 


geometrical and dynamical cosmic observations, we obtain 
log(|/ R o — 1|) = —5.71+Qjg which is consistent to the previ¬ 
ous result obtained by MGCAMB and EFTCAMB. 
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